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The luminosity function of diverse satellite galaxy systems 
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ABSTRACT 

The high-resolution, SPH galaxies of the McMaster Unbiased Galaxy Survey (MUGS) 
are used to examine the satellite systems of sixteen model host galaxies. Each galaxy 
has a different mass, angular momentum and merger history that yield a rich set of 
satellite luminosity functions. With new observations of distant satellite systems, we 
can compare these luminosity functions to satellite systems beyond the Local Group. 
We find that the luminosity functions of our simulations compare well to observations 
when the luminosity functions are scaled according to host mass. We use the recently- 
found relationship betwe en dwarf satellites and host mass in distant satellite systems 
(|Trentham fc T ullv 2009) t o normali ze a theoretical, complete luminosity function for 
the Milky Way (jKoposov et al.ll2008h . The luminosity function of satellites, expressed 
as a function of the host mass, is given by dN/dM v = 3.5M^ S \ x io 01M v-io.2^ w k ere 
mass is given in Mq . The mass of a host galaxy can be used to predict the number of 
dwarf satellites and even when considering spiral and elliptical hosts separately this 
relation holds. 
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■ 1 INTRODUCTION 
y—i [ 

^s^j Dwarf galaxies are the most common objects in the Uni- 
verse today l|Marzke fc da Co sta 1997). Their numbers pro- 
vide insight into many aspects of the evolution of the Uni- 
verse from the composition of dark matter to how star for- 
mation affects galaxy formation. Dwarf galaxies are flow- 
s' 1 ever difficult to observe because of their low surface bright- 
ness and a detailed census of dwarf galaxies on large scales 
remains a daunting challenge. Only for the most nearby 
systems is a complete counting of dwarf galaxies possible. 
Such accounting leads to the well-known discrepancy be- 
tween the lower observed number of satellites in the Local 
Group and the higher number predicted in nu merical simu- 
lations l|Moore et al.lll999l : iKlypin et al.lll999l ), also known 
as the missing satellites problem. 

Up until now, however, the observed data set has been 
limited. Two recent developments, observation of fainter 
dwarfs and the discovery of satellites outside the Local 
Group, for the first time yield data sets that allow for a 
study of the large-scale statistics of satellites galaxies. 

Large surveys, including the Sloan Digital Sky 
Survey (SPSS) (lAbazajian et all 120091) and PAn- 



dAS dMcConnachie et all 120081 iMartin et all 120091 : 



surface brightness levels and have revealed a new category 
of dwarf galaxies: the ultra faint dwarfs (Willma n et al.l 
120051 : iBelokurov et~aH 120071 : iKoposov et all I2008D . CuT 



rently these large surveys do not cover the entire sky and 
there may well b e ultra faint d w arfs t oo faint for current 
instrumentation. IKoposov et al. (2008) have calculated a 
theoretical and complete luminosity function of the Milky 
Way, based on observations of ultra faint dwarfs and it 
compensates for the incomplete sky coverage. 

Beyond ultra faint dwarfs within our own Local Group, 
dwarf satellites in other systems are being revealed as 
well. It is preferable to study satellites in varied environ- 
ments, as opposed to the Local Group exclusively, even 
though the Milky Way appears to be typical in terms of 
the a bundance of classical satellites (jStrigari fc Wechslerl 
2011). Observations of galaxy clusters have previously 
yielded the luminosity function s of cluster ga l axies (e.g. 
lAbel]||l977l ; IJones fc Jones1ll980l : iBmggeli et al.lll985l ) and 
more recently follow the mass limit down to resolve the 
satellites of the galaxies themselves. These include M101 
jHofmbergill950l ). the nearby Virgo Cluster (jSandage et al.l 
Il985l ). the Fornax Cluster dFerguson fc Sandagel 119881 ). 



IRichardson et alj 201ll ). bring observations to even lower 



the Ursa Major Cluster (|Trentham et all l200ll). NGC 



5846 l|Mahdavi et all 120051). NGC 1407 ( Trentham et ail 
120061 ). NGC 5371 l|Tullv fc Trenthaml 120081 ). M81 
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dChiboucas et all l2009h NGC1023 (ITrentham fc Tuuvl 
120091 ), th e Coma cluster lichiboucas et aLlbOllI ), the Antlia 
cluster (jSmith Castelli et alJ 120111 1. and other systems 
using the Sloan Digital Sky Survey l|Tollerud et alj l201ll ; 
lLares et al.ll201ll ). 

These two developments from the observational end 
have alleviated the traditional missing satellites problem by 
increasing the number of satellite galaxies. From the the- 
oretical end, simulations that take baryonic physics into 
account have also helped to close the the discrepancy be- 
tween the number of o bserved satellites and predicted dark 
matt er subhalos (e.g. iGovernato et al.l 120071: iMaccio et alJ 



20091; lOkamoto et all 120091 ; IWadepuhl fc Springell l20ld : 



Nickerson et al.ll201 lh ■ In these simulations, a large fraction 



of subhalos contain only trace quantities of baryons and ac- 
cordingly may only be observed through gravitational lens- 
ing and other gravity-based methods. However, other details 
still remain to be solved; for example, comparison of the 
faintest satellites is compromised by the limited resolution 
of the simulations, and there are discrepancies between the 
internal kinematics of observed satellites and the structure 
of the dark matter subhalo s they are predicted to inhabit 
l|Bovlan-Kolchin et alj|201ll ). 

When studying this new wealth of dwarf galaxies, both 
observed and simulated, it is also important to consider how 
the environment, chiefly the properties of the host galaxy, af- 
fect the satellites. Previous stu dies have focussed on brighter 
and more massive galaxie s (e.g. lChristlein|[2000l ; lNichol et al.l 
l2003l ; lBalogn" et alj |2004). while we are concerned with how 
environment affects the newly-observed satellite galaxies in 
dist ant systems. 

iTrentham fc Tuuvl {2009) studied the spiral-rich group 
of galaxies around NGC 1023, whose luminosity function 
has a characteristic magnitude —25.5 < M* < —24.0 and 
faint-end slope —1.22 < a < —1.14 (|Press fc Schechterl 
1 1974 ISchechtedll976T ). They combined their new findings 
with the host galaxies NGC 1407, NGC 5846, NGC 5353/4, 
M81, and the Local Group. The galaxies' masses range from 
a little over 10 12 MQto almost 10 14 Mq, and each contain 
15 to 250 satellites within the stated credibility limit of 
Mr ^ —11. They divided the satellites into giants and 
dwarfs at Mr — —17 and compared the total number of 
satellites in each category to host mass. For the dwarf satel- 
lites they found a tight correlation with host mass: 



AT S/r 0.9:L±0.11 

N d ex M host ; 

while giants had a looser correlation: 



(1) 



(2) 



This suggests that the number of dwarf satellites is more 
accurately predicted by host mass, while giant satellites are 
not plentiful enou gh to provide as rel iab le sta tistics. 

We combine ITrentham fc Tullvb d2009h d ata s et of 
dwarf mass to host mass with Ikoposov et all s (|2008T l the- 
oretically complete luminosity function for the Milky Way 
to establish a final luminosity function that scales with the 
mass of a host galaxy. We will verify that our simulations 
are well-described by this, and also find that the relation 
remains the same for spiral and elliptical hosts. 

In section ijHwe briefly review the MUGS simulations 
and halo finding method; in jj3]we introduce the sixteen host 
systems studied and detail the luminosity function of each; 



in jj4]we compare our satellite luminosity functions to the 
Press-Schecheter function and the Trentham-Tully relation 
to justify its use; and our conclusions are 33 



2 METHOD 

We analyse the subhalos of sixteen gala xies from the McMas - 
ter Unbiased Galaxy Survey (MUGS) l|Stinson et alj|20ich , 
a sample of M* galaxies simulated at high res olution. MUGS 
was run using the SPH code GASOLINE (|Wadslev et al.l 
|2004). GASOLINE includes low-temperature metal cooling 
( Shen et al]|2010l ) , UV background radiati on, star formatio n 
that models the Kennicutt-Schmidt Law (|Kennicuttl fl998 ) . 
and physicall y-motivated stellar f eedback from the "blast- 
wave model" l|Stinson et al.l [2006!) ■ The met al cooling grid 
is con structed using CLOUDY (version 07.02 iFerland et al.1 
(1998)), assuming ionisation equilibrium. A uniform ultravi- 
olet ionising bac kground, adopted f r om H aardt & Madau (in 
preparation; see lHaardt fc Madau (199(1)), is used in order 
to calculate the metal cooling rates self-consistently. The UV 
starts to have an effect at z ~ 9 .9. With these prescriptive 
elements, iNickerson et al.l I^Oll! 1 ) found that the resulting 
model galaxies do not exhibit the missing satellite problem. 
This results from a combination of early UV heating, ram 
pressure stripping, tidal stripping and stellar feedback that 
substantially modify the mass-to-light ratios of the model 
galaxies' satellites. This is a prediction of the simulations 
and does not involve any specific tuning in the galaxy model. 

We evolve several 50 h~ J Mpc volumes of a WMAP3 
ACDM universe (H =73 km s _1 Mpc " 1 , D. m =0.24, 



fi A =0.76, JW„on=0.04, and a 8 =0.79) (jSpergel et alj|2007h ■ 
From these galaxies we choose a random selection with halo 
masses between » 5 x lO n M and » 2 x 1O 12 M that 
did not evolve within 2.7 Mpc of a structure more massive 
than 5.0 x 10 11 Mq. The sample is unbiased with regards to 
angular momentum, merger history, and less massive neigh- 
bours and it is desired that the sample will reproduce the ob- 
served spread in galaxy properties. The selected galaxies are 
then re-simulated with the commonly-used zoom technique, 
which adds high resolution dark matter and baryons in the 
region of interest, while maintaining the periferal galaxies 
at low resolution to provide the appropriate tidal torques. 
The initial dark matter, gas and star particle masses are 
1.1 x 1O 6 M , 2.2 x lO 5 M and 6.3 x lO 4 M respectively. 
Each type of particle uses a constant gravitational softening 
length, 310 pc. A fu ll description of MUGS can be found in 
IStinson et all (|201Ch . 

In order to identify a host galaxy and its subhalos, we 
use the Amiga Halo Finder (AHF) (|Knollmann fc Knebd 
2009). AHF is based on the spherical overdensity method to 
identify density peaks using an adaptive mesh algorithm. 
AHF cuts out halos (and subhalos) of identified density 
peaks using isodensity contours. A simple unbinding pro- 
cedure is used to determine whether the particles are gravi- 
tationally bound to a halo or its subhalos. 



3 SIXTEEN SYSTEMS OF SATELLITES 

MUGS consists of sixteen galaxies simulated at high reso- 
lution each with identical baryonic physics. Their masses, 
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Galaxy 


Mass 
(10 n M Q ) 


Class 


(kpc) 


N gas 
(10 5 ) 


(10 6 ) 


(10 5 ) 


Mass sa t 
(10 10 M0) 


N sa t 


N, 

/ n ft'i 


g7124 


5.0 


E 


165 


1.4 


1.2 


3.7 


7.8 


37 


9 


g5664 


5.7 


S 


173 


2.0 


1.1 


4.3 


3.3 


43 


8 


g8893 


6.7 


E 


182 


2.2 


1.4 


5.0 


2.5 


57 


10 


gl536 


7.5 


S 


190 


2.7 


1.4 


5.7 


1.5 


56 


7 


g21647 


8.8 


s 


200 


3.0 


1.8 


6.7 


14 


63 


10 


g22795 


9.2 


E 


203 


3.1 


1.5 


7.1 


3.5 


76 


11 


g22437 


9.5 


E 


206 


3.8 


1.7 


7.2 


7.3 


72 


9 


g422 


11 


S 


218 


1.0 


2.3 


8.5 


32 


108 


16 


g3021 


11 


E 


218 


3.9 


2.5 


8.5 


7.8 


110 


17 


g24334 


12 


s 


221 


3.7 


2.5 


8.9 


18 


113 


26 


g28547 


13 


s 


226 


4.0 


2.9 


9.1 


26 


107 


18 


g4720 


13 


s 


229 


5.2 


2.1 


10 


17 


151 


27 


g25271 


14 


E 


233 


4.6 


2.4 


11 


3.3 


96 


12 


gl5784 


15 


s 


240 


5.3 


2.6 


12 


11 


96 


21 


g4145 


15 


s 


239 


5.7 


2.8 


1.1 


19 


130 


29 


gl5807 


23 


E 


276 


8.7 


4.0 


17 


11 


161 


33 



Table 1. The attributes of the sixteen MUGS host galaxies at redshift zero, arranged in order of increasing mass: mass, the classification 
(.Elliptical or Spiral), the virial radius as found by AHF, the number of gas, star and dark particles, the total mass of satellites found by 
AHF with fifty or more particles at redshift zero, and the number of satellites and the subset of those that are luminous. 



morphologies, virial radii and the number of gas, star, and 
dark particles at redshift zero are given in Table [3] illustrat- 
ing the diversity of MUGS. The total mass of satellites in 
the system found by AHF with fifty or more particles, the 
number of satellites and the number that are luminous at 
redshift zero are also included. g8893, g422, g3021, g28547, 
g4720, g4145, and gl5807 are presen ted here for the firs t 
tim(Q, while the rest were presented in lStinson et alJ l|2010l ). 

For reference, current estimates of the ma ss of 

the Milky Way are around 10 12 M g [|Klvpin et alJ 1200 
iBovy et al.ll2012l) and and for M31 are around 5.0 x 10 12 
M0 (|Foreman-Mackevll2010l ). Finding these masses is still a 
very active and ongoing field of research. 

Details of how we calc ulated the luminosity for each 
subhalo can be found in Nickerson et all (l201lh. using 
the initial mass functions of iKroupa et al.l (1 19931) and 
the luminosity gri d from CMP 2.1 l|Leitherer et all Il999l ; 
iMarieo et aT1l2008l ). 

Figure [1] shows the cumulative V-band luminosity func- 
tion of the subhalo populations of the sixtee n host galaxies, 
ordere d by mass, at z = 0. Also shown is the lTollerud et aTl 
(2008) data for the Milky Way that includes both the 
classical satellites and the new ultra-faint dwarf galaxies 
(which are fainter than the resolution of our simulations) 
and the luminosity func tion of M31. The classical M31 
satellites are taken from iMateol (jl998j}_ with the addition 
of newly-discovered ultra faint dw arfs ( Brasseur et aljfeoill ; 
ISlater et al.ll201ll ; iBell et alj|201ll ). 

It is believed that the set of ultra faint dwarfs observed 
in the Local Group is incomplete because of their low surface 
brightness and incomplete sky cover age of the Sloan Digi- 
tal Sky Survey. IKqposov et al.l (|2008l ) provides a theoretical 
function that would represent a complete set of subhalos 
with V-band magnitudes from -2 to -11 for the Milky Way 
given by 



1 Images of these galaxies can be found at 
http:/ /mugs.mcmaster.ca/ 



dN 
dM v 



= 10 x 10 



0.1(M v +5) 



(3) 



where N is the number of satellites that have a magnitude 
of My or brighter. However, in order to look at the lumi- 
nosities for systems of satellites in non-Milky Way-like hosts 
we need to readjust the normalization, i.e. the total number 
of satellites expected in different sized galaxies. For this w e 
use the Trentham-Tully relation (|Trentham fc Tullvl [20091 '): 



logA d = -10.2(±1.4) +0.91(±0.11)logM hos 



(4) 



where Nd is the total number of dwarf galaxies with ri- 
band magnitudes between -11 and -17 for host of mass 
Mhost in Mq. We only use the dwarf satellites for this and 
not the giants, for reasons that are detailed in jQ ] Some 
spec ulation is involved he re because iKoposov et all (|2008h 
and Tollerud et all (T2008I) have differing magnitude limits, 
iTollerud et al.l ( 20081 ) being limited by low surface bright- 
ness, which is why we will need to test this in §3] and ensure 
that it works. Because the range in the Trentham-Tully re- 
lation is given in the R-band, while the V-band is used for 
the Koposov function, we need to convert the range. From 
our sample of 280 luminous satellites in MUGS, we find that 
V — R — 0.51±0.01. The final expression for the luminosity 
function as a function of the host mass and V-band magni- 
tude is: 



3.5MS x io°- 1M v-io.2 



(5) 



dN 
dM v 

and this is also shown in Figure [T] alongside the luminosity 
functions of the MUGS simulations. 

The Milky Way and M31 luminosity functions have dif- 
ferent shapes. The Milky Way remains as a steady power 
law for high luminosity satellites. M31 has a knee, showing 
a higher number of high luminosity satellites. Correspond- 
ingly, a few of our galaxies also match the M31 shape bet- 
ter than the Milky Way's, notably g422, g3021 and gl5784. 
On the other hand, g8893, g22795, and g4145 match the 
shape of the Milky Way luminosity function, as given by 
Toll erud et al. (2008) closely. However, most other galaxies 
display an excess knee in high luminosity satellites compared 
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g7124 
5.0X10"M 




g5664 
5.7X10"M 




g8893 
6.7X10"M 




g1536 
7.5X10"M 




g21647 
8.8X10 1, M 



g22795 
9.2X10 1, M 





g22437 
9.5X10 1| M 



g422 




1.1X1 12 M 




g24334 
1.2X1O 1z M 



g28547 
1 .3X10 U M 





g4720 
1.3X10'^M 




g2527 
1.4X10"M 




g15784 



g4145 
1 .5X1 O'^M 




g 15807 
2.3X1 12 M f 



■10 -15 



■10 -15 



10 -15 



■10 



Figure 1. The cumulative V-band luminosity functions of the subhalo populations of the MUGS ho st galaxies at z = 0, in solid black, 
ordered by host mass, with the host mass given at the bottom. Observational data of the Milky Way llTollerud et al, |2008|) is in dashed- 
dotted red, M31 is in triple-dot-dashe d purple llMatec Il998l : iBrasseur et al.ll201ll : Islater et al.ll201ll iBell et al.ll20li) and our E quation 
[5](KTT) is in dashed blue taken from lKoposov et"HrtT 2009?) . normalized to our Equation Q] taken from iTrentham fc Tullvl |2009). . The 



ultra faint dwarfs as observed in the Milky Way go much fainter than the resolution limit of our simulation. 



to the Milky Way. We will now examine the total number 
itself. 



Of the galaxies that are less massive than the Milky 
Way and go up to its mass, ranging from ~ 50% to 95% of 
its mass, g7124 and g8893 have a similar cumulative number 
of satellites compared to the normalized power law Equation 
[5] g5664, g!536, g21647, g22795, and g22437 are relatively 



satellite-poor, though are still within an order of magnitude 
of what is predicted. 

The Milky Way-like-in-size galaxies and slightly more 
massive (g422, g3021, g24334, g28547, g41720, g25271) 
have very sim i lar cu mulative luminosity functions to the 
iTollerud et alj (|2008l ) function for their brightness limit, or 
are near to the Milky Way's function. The MUGS curves do 
not appear to dramatically change shape as they approach 
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the resolution limit, which suggests that the total number 
of low luminosity satellites is relatively robust. 

Of the heavier galaxies that exceed the Milky Way's 
mass, gl5784 is the only one of these that meets M31's satel- 
lite count and g4145 and gl5807 match Equation[5] Our high 
luminosity satellites are probably more luminous than is re- 
alistic, but they are within the resolved range and should 
still contribute to the total cumulative number. 

The importance of having a sample of galaxies in dif- 
ferent environments is demonstrated by the variance in lu- 
minosity functions across galaxies of similar mass. Overall, 
in terms of cumulative number, and within the resolution 
limit, our galaxies do not suffer from the order of magnitude 
missing satellites problem as is evident from Equation [5] In 
spite of the differing environments, a trend still holds to pre- 
dict the number of satellites based purely on the host's mass. 
Only the satellites below the Milky Way's mass have slightly 
fewer satellites than expected, while every Milky Way-mass 
and more massive host meet this trend. This might be an 
effect of small scale statistics for the lower mass hosts and 
with higher resolution we expect they too should follow the 
trend closely. 



4 THE ROBUSTNESS OF DWARF 

SATELLITES OVER GIANT SATELLITES 

Figure [2] shows the cumulative luminosity function for all 
the satellites of the sixteen MUG S galaxies. It includes a fi t 
to the Press-Schechter function l|Press fc Schechterl Il974l ). 
adjusted to an intercept such that the brightest s atellite has 
a cum ulative number of 1. In cumulative form (Schechter 
119761 ) the function is, 

N a = T{a + l,L/L*)n* (6) 

where their N e is the number of galaxies expected to have 
a higher luminosity than L, a is the faint-end slope, L* is 
the characteristic luminosity at which the faint end begins, 
n* is the richness parameter, and F is the incomplete upper 
gamma function. 

Most of our galaxies display an excess of high- 
luminosity satellites, making the fit to the Press-Schechter 
function somewhat problematic. We were nevertheless able 
to obtain a fit that brackets the high-luminosity bend: 
a + 1 = -5.48 ± 0.67 x 10" 2 and M* = -21.3 ± 0.2. It is 
already clear from Figure [2] that our dwarf satellites fit the 
Press-Schecter functional form much better than the giant 
satellites. 

Figure [3] shows the number of lumi nous satellites ver- 
sus th e mass of the host galaxy. As with lTrentham fc Tullvl 
(2009) we do not consider satellite galaxies fainter than 
Mr = —11, and we split the dwarfs from the giant satel- 
lites at Mr — —17. The power law fits for this relation are 
shown in Figure [3] We obtain for dwarfs, 

iV d oc Mt 2 ^ 2 (7) 
and for giants, 

N a txMZ? A . (8) 

Within the errors, our trends for the dwarfs and giants fall 
within the Trentham-Tully relations given in Equations [T] 
and [2] We find, just as they do, that the giant satellites 



VII 




-6 -8 -10 -12 -14 -16 -18 -20 
Mr 

Figure 2. The cumulative R-band luminosity function for all 
satellites of the MUGS hosts in solid blue. The Press-Schechter 
functional fit is in dashed red. 

have more scatter when fit to a power law than the dwarf 
satellites. 

We investigated splitting the hosts by morphology and 
exploring the difference in the power law between spiral and 
elliptical hosts in Figures [4] and [5] While the power law 
for the number of dwarf satellites (1.3 ± 0.4 for spirals and 
1.2 ± 0.2 for ellipticals) remains nearly the same, the power 
for the giant satellites changes drastically between the two 
morphologies (1.8±0.5 for spirals and 0.42±0.58 ellipticals). 
This suggests further that the small-number statistics of gi- 
ants are more affected by environment, while the number of 
dwarf satellites is robust against morphology type. 

Dwarf satellites are more useful for finding a universal 
scaling relation for the luminosity function of satellites ac- 
cording to host mass as in Equation [5] These observations 
verify the trend shown in our simulations, in that the mass of 
the host predicts the number of dwarf satellites, and hence 
within resolution our simulation does not suffer from the 
missing satellites problem. This is an important loadstone 
for studying luminosity functions and satellite systems out- 
side of the Local Group. 



5 CONCLUSIONS 

We compare the satellite populations in simulations of six- 
teen high resolution galaxies with observations to see how 
satellite populations vary as a function of host mass. Our 
sample of hosts contains a wide range of masses and mor- 
phologies and are therefore interesting to compare to the 
new studies and statistics of satellite galaxies in diverse sys- 
tems beyond the Local Group. Just as there are observations 
of diverse galaxies, there needs to be diverse simulations run 
with the same physics that contain more than Milky Way- 
like galaxies. 

The luminosity functions of our galaxies scale well with 
observations. It shows that with baryon physics, our simu- 
lations no longer suffer from the order of magnitude miss- 
ing satellites problem. For host galaxies as massive as the 
Milky Way and more massive, the scaling fits almost ex- 
actly. Though our hosts with masses lower than the Milky 
Way do not quite have as many satellites as predicted, they 
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10 - 



■ giants (M R ^- 1 7) 

• dwarfs (-17< M„5-11) 



N H 



N g!ont K M hos. 



1.0 ±0.4 
1.2 ±0.2 



10" 

Host Mass (M ) 

Figure 3. The number of satellites versus host mass. The satel- 
lites are divided at R-band magnitude -17 between giants (red 
squares) and dwarfs (blue circles). The power law fit for the gi- 
ants is in dotted-dashed red, and the fit for the dwarfs is in dashed 
blue. The giants are offset to the right slightly so that if they over- 
lap in number with the dwarfs for a particular galaxy both points 
can still be seen. 



10 - 



■ elliptical host, giants 
• elliptical host, dwarfs 

M ~ u 0.42±C 

N„,„„. « M h „ e , 



N H 



10" 

Host Mass (M ) 

Figure 5. The number of satellites versus host mass for elliptical 
host galaxies only. The satellites are divided at R-band magnitude 
-17 between giants (red squares) and dwarfs (blue circles). The 
power law fit for the giants is in dotted-dashed red, and the fit 
for the dwarfs is in dashed blue. The giants are offset to the right 
slightly so that if they overlap in number with the dwarfs for a 
particular galaxy both points can still be seen. 



10 



■ spiral host, giants 
• spiral host, dwarfs 



10" 

Host Mass (M ) 

Figure 4. The number of satellites versus host mass for spiral 
host galaxies only. The satellites are divided at R-band magnitude 
-17 between giants (red squares) and dwarfs (blue circles). The 
power law fit for the giants is in dotted-dashed red, and the fit 
for the dwarfs is in dashed blue. The giants are offset to the right 
slightly so that if they overlap in number with the dwarfs for a 
particular galaxy both points can still be seen. 



were still only short by three or four satellites, and probably 
with higher resolution would suffer less from small number 

stat istics. 

iTrentham fc Tullvl 's (|2009l ) study concluded that dwarf 
galaxy populations (—17 < Mr < —11) correlate to their 
host halo's mass, while giants (Mr < —17) do not. They 
explored the space of host masses from 10 12 to 10 14 Mq 
and found the relation Ndwarf oc M^f^ ' 11 , where mass is 
given in Mq. We find that that the Trentham-Tully rela- 
tion describes our simulated satellite systems even when we 
split our hosts between spirals and ellipticals. The number 
of dwarf satellites a host has is largely dependent on the host 



mass and does not change with morphology type, whereas 
the number of giant satellites might be affected by many 
factors. There were not enough giants on which to found 
solid statistics. 

When all satellite galaxies from all sixteen MUGS sim- 
ulations are stacked together, their luminosity function fol- 
lows a Press-Schechter function with faint-end slope a + 1 = 
—5.48 ± 0.67 x 10 -2 and characteristic luminosity M* = 
-21 3 ± 0.2. 

IKopqsqv et al.l l|2008h provided a luminosity function 
for the Milky Way that takes into account the incomplete 
sky coverage of large surveys and dwarfs too faint for detec- 
tion. We use the knowledge of dwarf satellites' robustness 
to adjust the Milky Way's function in order to incorporate 
dependence on host mass, as in the Trentham-Tully rela- 
tion for dwarfs: dN/dM v = 3.5iW£ 9 s < x W°- 1M v - 10 - 2 , where 
mass is given in Mq. The MUGS hosts compare much more 
favourably to this new, scaled luminosity function than they 
do to the Milky Way's or M31's luminosity functions alone. 
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